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Statistical properties of exoplanets. 

I. The period distribution: constraints for the migration scenario 
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Abstract. Interesting emerging observational properties of the period-mass distribution of extra-solar planets are 
discussed. New recent detections confirm the already emphasized lack of massive planets (mj sini > 2Mj up ) on 
short-period orbits (P < 100 days). Furthermore, we point out i) a shortage of planets in the 10-100 day period 
range as well as ii) a lack of light planets (7712 sin i < 0.75 Mj up ) on orbits with periods larger than ~ 100 days. The 
latter feature is shown not to be due to small-number statistics with Monte-Carlo simulations. These observational 
period-related characteristics are discussed in the context of the migration process of exoplanets. They are found 
to be in agreement with recent simulations of planet interactions with viscous disks. The observed valley at a 
few tens of days in the period distribution is interpreted as a transition region between two categories of planets 
that suffered different migration scenarios. The lack of light planets on longer-period orbits and the corresponding 
intriguing sharp limit in mass is tentatively explained by the runaway migration process recently studied by 
Masset & Papaloizou (2003) . The observed properties also have implications for the observation strategies of the 
on-going surveys and of future higher-precision searches. 
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1. Introduction 

The most remarkable feature of the sample of known 
extra-solar planets is undoubtedly the variety of their 
orbital characteristics, which challenges the conventional 
views of planetary formation. Amongst the most pecu- 
liar candidates are the giant planets orbiting very close 
to their parent stars , in contrast to the prediction of the 
standard model (e.g. lPollack et aiT l996) that they formed 
first from ice grains 1 in the outer region of the system 
where the temperature of the stellar nebula is not too high. 
The accommodation of this scenario to the present obser- 
vations requires that the planets undergo a subsequent 
migrati on process bringing them cl ose to the central star 
(see e.g. lLin et al.lll99fjlWardll997h. Alternative p oints of 
view invoke in-situ formation l)Wuchterlll2000albr) , poss i- 
bly triggered through disk instabilities l)Bossll200ll 12002^ . 
Note however that, even in such cases, subsequent disk- 
planet interactions leading to migration are expected to 
take place as soon as the planet is formed. 
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1 Such grain growth provides the supposed requisite solid 
core around which gas could rapidly accrete, over the lifetime 
of the protoplanetary disk (~ 10 7 y) 



The number of known extra-solar planets exceeding 
100, a statistically significant sample is now available 
from which we can deter mine meaningful distributions of 
plane tary characteristics l|lldrv et al.1 2003at iMarcv et alJ 
2003), and so try to point out useful constraints for the 
planet-building models and then possibly discriminate be- 
tween the different proposed scenarios. 

In a new series of papers, we will try to emphasize the 
emerging properties of planet-host stars and characteris- 
tics of the different orbital-element distributions of exo- 
planetary systems and discuss their implications for our 
understanding of planetary formation and evolution. This 
paper (Paper I) will be dedicated to the period/separation 
distribution of exoplanets. Sharp features in the period 
versus mass diagram start to emerge and provide strong 
observational constraints for the migration scenario. More 
specifically, we re-discuss the lack of m assive planets on 
short-period orbits recently pointed out llUdrv et alJ l2 002; 
IZucker fc Mazehl l2002t IPatzold fc Rauerll2002|) .~We also 
emphasize the clear emergence of a planet shortage in 
the 10-100 d period range (a period valley) and we point 
out a sharp lower mass limit for the detected planets on 
"longer" -period orbits (P > 100 d). The later part of the 
paper will be dedicated to a discussion of the proposed ex- 
planations of these important observational findings and 
the possible implications for the future surveys. 
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Fig. 1. Upper panels: Minimum masse s versus periods for k nown exoplanet candidates. In the left panel, filled squares 
indicate planets in binaries fPaperlll: llldrv et alJEof)3flch whereas circles are used for planets around single stars. 
The point in " Q" represents HP 162020 for which synchronisation arguments indicate a probable mass in the brown- 
dwarf regime llUdrv et al.ll2002t) . In the right panel, only planets orbiting single dwarf stars are represented (i.e. 
planets in binaries and planets orbiting evolved stars have been discarded). A different coding is used for massive 
(m,2 sini > 2 Mj up ; filled symbols), intermediate-mass (Tasini between 0.75 and 2Mj up ; open circles), and lighter 
(7712 sini < 0.75 Mj up ; open triangles) candidates. The components of the possible multi brown-dwarf system HD 168443 
are linked by a dotted line. The dashed and dotted lines in the panels indicate limits at P = 100d (vertical), at 
7Ti2 sin i — 2 Mj up (horizontal left), or at mi sini = 0.75 Mj up (horizontal right). Lower panels: Period distributions 
of the planets orbiting single dwarf stars (shown in the upper-right panel) for different mass regimes: light-shaded 
histograms are for the lightest planets [m-i sini < Md OW n = 0.75 Mj up ), open histograms are for intermediate masses 
(between Md own and M up - right: M up = 2Mj up ; left: M up = 4Mj up ), whereas the dark-shaded histograms are for the 
more massive candidates (7712 sini > M up ). 



A second paper in this series ((Santos et al.1 [2003, 
Paper II) deals with the metallicity of the stars harbour- 
ing planets, providing accurate spectroscopic parameters 
for the most recent candidates, confirming known global 
properties and a lso pointing out new eme rging features. 
In a third naner. lEggenberger et alJ l|2003l Paper III) ex- 
amine the effect of binarity on planet formation and evo- 
lution. The presence of a close stellar companion is found 
to bring more massive plan ets close-in, in agreement with 
simulations bv lKlevI l)200l|) . 



2. Observational emerging properties 

2.1. No massive planets on short-period orbits 

In a recent discussion of the st atistical propertie s of mas- 
sive planets versus lighter ones ijlTdrv et al.ll2002ft . we em- 
phasized a shortage of massive planets on short-period 
orbits, based on the sample of about 80 exoplanets known 
at that time. This feature was al so simultaneously pointed 
out by IZucker fc Mazehl §002). These authors further- 
more verified its statistical significance and very interest- 
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ingly examined the possible influence of binarity on the 
mass-period relation of exoplanets, an indication of poten- 
tial different formation and evolution processes for planets 
in binaries and planets around single stars. 

With now more than 100 very low-mass candidates 
(with m2sini <18Mj up ), the lack of short-period mas- 
sive planetary companions becomes even clearer, as seen 
in Fig.Q] When we neglect the multiple-star systems 2 
(upper-right panel) in which the planet ary formation or 
evolution could follow di fferent paths l|Zucker fc Mazebl 
l2002t lUdrv et al.ll2003blcl Paper III), a complete void of 
candidates is observed in the diagram for masses larger 
than ^2Mj up and periods smaller than ~ 100 days 3 . The 
only remaining point is HD 168443 b, member of a possible 
multi brown-dwarf system I Marcv et al.ll200lt llldrv et afl 
120021) . 

Several processes have been proposed to explain the 
lack of massive planets on short-period orbits. In the con- 
text of the migration scenario, they mainly follow two dif- 
ferent approaches: i) type II migration (after a gap opens 
in the disk) is shown to be less effective for massive plan- 
ets i.e. massive planets stay farther out than lighter ones, 
or ii) when the planet reaches the central regions, some 
process related to planet-star interactions provokes mass 
transfer from the planet to the star - decreasing the mass 
of the former - or leads massive planets to fall into the 
central star. These approaches will be discussed in more 
detail in Sect. El 

The possibility that multi-planet chaotic interactions 
send the lighter candidates in the inner regions (or out) 
of the system whereas the massi ve ones stay i n the 
outer part may also be invok ed ijRasio k. Fordl Il996t 
IWeidenschilling k. Marza"rilll996j) . The frequency of plan- 
ets endin g very close to th e central star seems however to 
be small llForr] et ahlEofil. 



2.2. A period valley between 10 and 100-d periods 

In the lower panels of Fig.^ and upper panel of Fig-El 
another very interesting feature is coming out of the pe- 
riod/separation distribution. We observe a shortage of 
planets with periods between roughly 10 and 100 days. 
This observational property is mainly due to the light 
candidates (rri2 sin i < 2 Mj up ) . As seen in the previous sec- 
tion, massive planets orbiting single stars are almost exclu- 
sively found on longer-period orbits. They probably form 
and stay further out (see discussion in Sect. 03). On the 
other hand, lighter planets are found at all distances from 
their star 4 . This mass-dependent behaviour is also clearly 




2 As determined by dedicated adaptive optics programmes 
or unveil ed by previous s pectroscopic measurements (see 
Paper III; lUdrv et alJl2003bl lcl for a review) 

3 The 4 known candidates of planets orbiting an evolved star 
have also been discarded. They however have no influence on 
the results of this study 

4 except the lightest ones (rri2 sini < ~ 0.75 Mj up ) at periods 
longer than ~ 100 days (see next section) 
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Fig. 2. Upper panel: Distribution of separations of the 
known exoplanets orbiting single dwarf st ars (log scale) 
compa red to the similar distribution of the iTrilling et all 
(2002) simulation end-states (dashed line). The planets 
more massive than 2Mj up are indicated by the shaded 
part of the histogram. The planet shortage in the 0.06 < 
a < 0.6 AU separation range (10-100 d) comes out very 
clearly. Lower panel: Cumulative functions of log a for the 
2 mass regimes (limit at 2Mj up ). The preponderance of 
light planets close-in and the location of massive plan- 
ets further out is clearly emphasized. The Kolmogorov- 
Smirnov probability for the 2 distributions to come from 
the same population is 6.7 ■ 10~ 4 . 



illustrated by the cumulative functions presented in the 
lower panel of Fig.[3for two planetary mass regimes. 

The features appearing in the planetary period distri- 
bution - peak at short periods and rise at intermediate 
periods - seem significant and not due to observational 
biases. The peak at short period is formed by the pile- 
up of migrating planets, stopped close to the central star. 
As seen in Fig.Q] it is almost exclusively composed of the 
lowest mass planets [m^ sin i < 0.75 Mj up ) and not due to 
observational biases as the detection limit at 10ms _1 is at 
P~4000 days for m 2 sin i = 0.75 M Jup (with M\ — 1 M 
and e — 0. The effect of Jupiter on the Sun is 12ms _1 ). 
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The rise of the distribution at longer periods is now emerg- 
ing thanks to the increase of the timebase of the radial- 
velocity surveys. More and more long-period planets are 
being detected. This rise is significant in the sense that, for 
a given mass range, similar planets on shorter-period or- 
bits would have been easier to detect than the actual ones. 
The observed valley in the period distribution can then 
just be seen as a transition region between two categories 
of planets which suffered different migration behaviours. 

2.3. A sharp mass transition in the migration process? 

A third important feature emerging from the mass- 
period diagram is the apparent lack of very light planets 
(to2 sin i < 0.75 Mj up ) with periods larger than ^100 days. 
This feature is already visible in the upper-right diagram 
of Fig.Q (limited by the dotted line). It becomes obvious 
when the planetary minimum mass is displayed with a log 
scale (Fig.0. 

Of course this feature could be related to the obser- 
vational bias inherent to the radial-velocity technique for 
planet search that makes the detection more difficult for 
distant and/or lighter planets. To obtain a quantitative 
idea of the effect of this bias, we have also plotted in 
Fig-El lines indicating the locations of the radial-velocity 
signals (semi-amplitude K of 30, 10 and 3 ms -1 ), expected 
on a solar-mass star due to planets on circular orbits 
with given minimum masses and semi- major axes 5 . The 
line for K = 10 ms" 1 gives the approximate 3er-limit of 
the present most precise surveys reaching precis ions of 
-3 ms" 1 (e.g. lVogt et"aIll2000l:lQueloz et aljbdoijl . Some 
known candidates have been found very close to this limit, 
at low mass and small distance or at high mass and 
large distance 6 . However, we clear see an "intermediate" 
region (hatched area), delimited by m-i sin i = 0.75 Mj up , 
P=100d and K = 10 ms -1 , where no planetary candi- 
dates have been found. A striking feature of this area is 
its very sharp limit in mass. 

Is this empty region just due to the mentioned obser- 
vational bias and the small-number statistics? We have in 
particular to worry about the observed effect coming from 
the shortage of candidates with periods in the 10-100 day 
range, as pointed out in the previous section. So, to check 
the statistical significance of the emptiness of the hatched 
area in the figure (i.e. to see whether it mainly comes from 
the small number of known planets), we performed Monte- 
Carlo simulations. The idea is to obtain 10 6 realizations 
of the diagram, reproducing the actual observational bi- 
ases of the known sample, and estimate from them the 
probability to have no point in the considered zone. 



5 The primary mass and eccentricity effects are not expected 
to be important as the radial-velocity semi-amplitude scales 
with (1 - e 2 )-^ 2 M- 2/3 

6 We can also note that these candidates "close" to the de- 
tection limit do not present special values of the eccentricity 
and primary or planetary masses 
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Fig. 3. Mass-separation diagram for the known exoplanet 
candidates. The dotted lines illustrate the radial-velocity 
semi-amplitude (3, 10 and 30 ms -1 ) expected on a solar- 
mass star due to planets on circular orbits with given min- 
imum masses and separations. The shaded area empty of 
planets is shown not to be due to small number statis- 
tics (see text). Planets in binaries are indicated by open 
symbols. 



There are basically two ways of producing samples re- 
sembling to the observations: either i) take realistic plane- 
tary parameter distributions and model the observational 
biases or ii) draw the test samples from the actually ob- 
served distributions in which the acting biases are built 
in. As the "real" distributions are not known a priori, we 
have adopted the second approach: 

i) We randomly select the orbital eccentricities, the pri- 
mary masses, the planet minimum masses and separations 
from the actual observed distributions of e, Mi, TTi^sini 
and log a of the detected exoplanets. This ensures in par- 
ticular that we are taking into account the effect of the 
orbital eccentricity and of the stellar mass on planet de- 
tection. It also produces distributions that reproduce the 
observed brown-dwarf desert and the period valley. 

ii) Since the distributions of separations are different 
for light and massive planets (Fig.0 , we do the selection 
for two mass regimes independently (limit at 4Mj up ). 

iii) We only use the distributions for planets orbiting 
single dwarf stars as 1 ) planets in binaries m ay have dif- 
ferent characteristics llZucker fc Mazer]l2002t Paper III) 
and 2) the post main-sequence stellar evolution may dras- 
tically change the fate of short-period planets. 

iv) We reject the events giving a radial-velocity 
semi-amplitude smaller than K\i m — 10 ms -1 . They would 
hardly be detected by the present surveys. A less restric- 
tive limit a 15 ms -1 is also considered to check the sensi- 
tivity of the results on the detection limit. 
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Fig. 4. Result on the statistical significance of the 
hatched zone in Fig.|3| based on 10 6 realizations of the di- 
agram with Monte-Carlo simulations. The figure presents 
the normalized histogram of the distribution of the num- 
ber of points, n, found in the test region i.e. the occurrence 
probability of n. The arrow location at n = indicates the 
position of the actual observational result. Its probability 
of occurrence only by chance is P(Q) — 0.00028. 

v) Each realization consists of a selection of 86 pseudo- 
planets, corresponding to the actual observations (64 
"light" and 22 "massive" planets orbiting single dwarf 
stars, as on February 2003). For each of them we count 
the number n of points in the "test zone" . Note that the 
detection limit at Kn m is no more as well defined as in 
Fig. El but it now takes into account the orbital eccentric- 
ity and primary-mass values. 

The normalized histogram of the distribution of the 
number of points n in the test zone from our 10 6 re- 
alizations is shown in Fig.0J The probability of having 
no detection is only 0.00028. A bootstrap of the pro- 
cedure with 100 times 10 5 realizations gives a value of 
0.00025 ±0.00006 (Fig.©, in agreement with the previ- 
ous result. Relaxing the detection limit to 15 ms^ 1 yields 
a probability of 0.004. It shows the low sensitivity of the re- 
sult on the chosen radial- velocity detection limit. We thus 
can conclude that the considered area is indeed empty of 
planets with a confidence level of about 99.97%. 

This conclusion is only valid if the observational bi- 
ases are correctly taken into account. Biases related to 
planetary parameter distributions and detection technique 
are by construction included in the simulations. However, 
other aspects play a role. For example activity-induced 
radial- velocity jitter may screen planet detection. Such an 
effect can be neglected if we suppose that the r^sinz, 
e, Mi and a parameters are not correlated with activity. 
Another concern relates to the difficulty encountered when 
trying to actually derive the values for the orbital ele- 



Fig. 5. Distribution of the P(0) probabilities of 100 boot- 
strap simulations (but with 10 5 realizations each) con- 
ducted to estimate the uncertainty on the value obtained 
in Fig.H They yields a value 0.00025 ± 0.00006 for P(0). 

ments, what is different from just detecting radial- velocity 
variability. Observation timing and phase coverage are 
then also important parameters. These remarks emphasize 
the caution needed when considering the above described 
results. We should see them more like a trend emerging 
from the data rather than a real proof. 

In the next section, we will propose a tentative expla- 
nation for this feature, ba sed on recent results o n run - 
away migration studied bv iMasset fc PapaloizoiJ l)2003h . 
We will also discuss in the last section of this paper the 
implications on the radial- velocity planet-search survey of 
the paucity of light planets with intermediate periods. 

3. Checking the theoretical approaches 

The different features pointed out above in the various 
presentations of the period/separation-mass diagram will 
allow us to check the processes proposed to explain the 
observed distribution of periods for extra-solar planet can- 
didates. 

3.1. Low migration rate for massive planets 

A first approach to explain the lack of massive planets 
on short-period orbits invokes a low migration efficiency 
for higher mass planets. This view is supported by recent 
simulations of single-planet migration. When the mass of 
the planet becomes of the order of a significant fraction 
of the characteristic disc mass with which it interacts, the 
inertia of the planet bec omes important and slows down 
the o rbital evolution (e.g. lTrilling et alJl998tlNelson et all 
2000). It also takes more time to open a gap in the disk 
for massive planets - eventually longer than the lifetime 
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of th e disk - and initiate type II migration ^Trilling et alJ 
12002ft . 

Using a model assuming a simple impulse approxima- 
tion for the type II migration, no ad-hoc stopping mech- 
anism in the cente r and neglecting type I migration 7 , 
Tril ling et al ] ll2002h show that, for identical initial disk 
conditions, close-in surviving pseudo-planets have smaller 
masses and distant surviving pseudo-planets have larger 
masses. They also find that planets forming farther out 
migrate less rapidly and that higher-mass planets migrate 
on longer timescales. Thus, there should be more massive 
planets at intermediate and large semi-major axes than 
found close-in. On the contrary, the population of short- 
period planets should be dominated by smaller-mass plan- 
ets, what is indeed observed: 

- The upper panel in Fig.[3] presents the logarithmic 
distribution o f exoplanet separations, superimposed to the 
iTrilling et all l|2002l) results. As already pointed out by 
these authors, if we take into account the observational 
bias penalizing long-period planets in the radial-velocity 
surveys, the theoretical and observational results agree 
fairly well. The observed peak at small separations is not 
reproduced by the simulations of Trilling et al. because 
they have not introduced in their study any ad-hoc pro- 
cess to stop the migration close to the central star. 

- Furthermore, the lower panel in Fig.Elgiving the cu- 
mulative functions for 2 mass regimes (limit at 2Mj up ) 
clearly shows that low-mass planets form the predominant 
population of short-period orbits whereas more massive 
planets are located further out. 

- Moreover, the lower panels of Fig.^ indicates that 
you have to go further and further out to find more and 
more massive planets. The hatched histogram represent- 
ing planets more massive than 4 Mj up rises up to 1000 days 
(left panel) - although such massive planets should be 
easier to detect at smaller periods - whereas the maxi- 
mum of the rise of the distribution of lighter-mass planets 
is around 400 days. If we decrease down to 2Mj up the 
limit between the 2 mass regimes, the mentioned features 
almost overlap (right panel), indicating a trend for the 
higher masses to be found at larger distances from the 
central star. This higher mass trend with period is ex- 
plicited in Fig.0 which displays the mean mass (filled 
circle) or higher mass (average on the 3 highest values; 
open circles) of planets in period smoothing windows of 
width log P[days] =0.2. The effect becomes visible around 
P = 40d and turns to be very effective around P— 100 d. 
Detection limits at 10 and 30 ms -1 are represented by the 
dotted lines. They show that observational biases related 
to the radial- velocity technique are not responsible for the 



7 The extremely short timescale of the type I migration rep- 
resents a serious theoretical problem for planet formation. 
However, suggestions have been recently made that co-orbital 
corotation t orque in thin visco us disks may decrease its dra- 
matic effect iArtvmowic and even stop or revert it for 
a large set of conditions ( Massed 120021) . This has however to 
be further examined 
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Fig. 6. Mean mass (filled circle) or highest mass (aver- 
age on the 3 highest values; open circles) of planets in 
period smoothing windows of width log P [days] = 0.2. A 
clear increase of planet maximum masses with period is 
observed, despite the fact that massive planets observed 
at a given period should be easier to detect at smaller 
periods. Detection limits at 10 and 30 ins -1 (Mi = 1Mq, 
e = 0) are represented by the dotted lines. 



observed feature. This effect is not surprising considering 
first that massive planets are preferentially formed in the 
outer regions where a larger amount of building material 
is available and second that higher-mass objects migrate 
on longer timescales. Furthermore, as the disk dissipates, 
the migration of more massive planets slows down earlier 
since their larger angular momentu m produces greater re - 
sistance to disk-induced migration l|Armitage et all2002|) . 



3.2. Runaway migration and lack of light planets with 
longer periods 

A recent study by iMasset fc Papaloizoul l|2003^ brings 
some insights for our understanding of the observed lack of 
small-mass planets (m? sin i < 0.75 Mj up ) on intermediate- 
period orbits (P > 100 d). These authors studied the effect 
of co-orbital corotation torque on migrating protoplancts 
in the case of massive protoplanetary disks for which the 
planet Hill radius and the disk thickness have comparable 
orders of magnitude. In particular they show that, if the 
mass deficit created by the radial drift with the planet of 
the material trapped in the co-orbital region is larger than 
the planet mass, the migration rate undergoes a runaway 
which can rapidly vary the protoplanet semi-major axis 
by a large amount (50% over a few tens of orbits). The 
authors state that: "this can happen only if the planet 
mass is sufficient to create a gap in its surrounding region 
and if the surrounding disk mass is larger than the planet 
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mass. This typically corresponds to planet masses in the 
sub-Saturnian to Jovian mass range embedded in massive 
protoplanetary disks". 

Mo reover, the simulations by iMasset fc PapaloizovJ 
l)2003(l show that the limit in planetary mass for the 
appearance of the runaway regime is very steep (their 
Figs. 12-14). The threshold conditions for runaway may 
vary widely depending on the disk mass, viscosity or thick- 
ness. It however depends relatively weakly upon the planet 
mass on the high-mass side i.e. their is a critical mass 
relatively well-constrained (-M cr i t ~ 1 Mj up ) under which 
runaway is likely provided that the protoplanetary disk 
is not too light weight, and above which runaway is im- 
possible (Masset, priv. comm.). As our observed limit at 
777,2 sin z ~ 0.75 Mj up is compatible with this critical mass, 
there could well be a link between the statistical observed 
property and runaway migration. 

In this picture, planets which happen to be massive 
enough (M>M cr it) before reaching the central regions 
would end up as the scattered population on the right 
of Fig.El whereas planets that never grow enough would 
remain very mobile and would eventually be flushed into 
the very central region. 

Note that another limit exist for runaway migration on 
the planet low-mass side. Thus, very light planets (typi- 
cally much lighter than Saturn) do not encounter such 
type of fast migration. That can be the case for the outer 
icy planets of the Solar System. 

3.3. Observed constraints for migration: summary 

In summary, putting everything together, we can try to 
draw a global picture for the planet migration behaviour 
in function of the planet mass, in single-star systems: 

i) The more massive planets (typically, masses larger 
than ~4Mj up ) form preferentially in the outer regions - 
where there is a large-enough material reservoir - and do 
not migrate much. None is observed within ~ 0.5 AU from 
the central star. 

ii) Intermediate-mass objects migrate more easily 
whatever the distance where they form, the migration rate 
depending on the local conditions (planet mass, disk mass, 
viscosity, etc). They are observed at all distances. 

iii) The lighter planets (masses from sub-Saturnian to 
Jovian) migrate easily, possibly undergoing runaway mi- 
gration bringing them in the central tenth of an AU from 
the star if the disk is massive enough. Most of them are ac- 
tually found in the very central regions. None is observed 
with P>~ 100 days although radial- velocity surveys are 
now precise enough to detect them. 

iv) We know, however, from our own Solar System that 
this is not true for the planets that are far away or much 
lighter than the known exoplanets. Thus, we can speculate 
that for the very light planets the migration efficiency is 
decreasing with distance or mass. This again fits into the 
the runaway migration scenario, the simulations showing 
the existence of a transition limit in mass, on the low- 



mass side, for the runaway to occur. Unfortunately, the 
precision achieved by the present surveys does not allow 
us to determine the separation or mass limits at which the 
migration slows down significantly. Good hope to answer 
this question is however brought by the future surveys at 
higher precision (like HARPS at 1 ms^. lPepe et al.l l2002) 
that will increase our detection sensibility by a large factor 
(curve at if = 3ms _1 in Fig.[3J). 

This observational view of planet migration may be 
changed by the influence of a perturbing stellar compan- 
ion (Paper III) that will make the massive planets migrate 
inwards faste r or make c loser-in intermediate-mass planets 
grow bigger llKlevll200lh . 

Finally, we have to stress that although theory suggests 
that it is possible to form massi ve planets closer to the 
stars than the ice-boundary limi t llPap aloizou fc Tera ueml 
Il999t iBodenheimer et all I2OO0I ISasselov fc Lecarl l200o|T 
within the migration scenario, there is no observational re- 
quirement for planets to form at separations smaller than 
a few AU. 

3.4. Central star-planet interaction 

From the central star - planet interaction point of view, 
several scenarios were proposed to explain the lack of mas- 
sive planets on short-period orbits as alternatives to the 
above developed arguments in the context of the migra- 
tion scenario. 

IPatzold fc Rauerl l|2002F) proposed to explain the lack 
of massive close-in planets by tidal interactions between 
the planet and its central star. If the planet orbital pe- 
riod is smaller than the stellar rotation period, tidal fric- 
tion will spin up the star whereas the semi-major axis of 
the planetary orbit is decreased, eventually reaching the 
Roche zone of the central star on a short timescale. These 
authors pointed out that massive planets spiral in much 
fas ter than other s , exp laining thus the observed trend. 

iTrilling et alJ l|l99Sj) had already addressed this prob- 
lem in also considering mass transfer between the planet 
and the star. Some of the pseudo-planets in their simula- 
tions reach the Roche lobe limit, start to lose mass and 
thus experience an outward torque from the Roche lobe 
overflow. They estimate the maximum mass of surviving 
planets to be around 2Mj up , corresponding to what is in- 
deed observed. 

Another approach could relate to the evaporation of 
the planet atmosphere when the planet comes close to the 
central star, decreasing thus its mass. Up to now, stud- 
ies based on the T e g of the heated close-in planets have 
shown that hydrogen evaporation was not efficient in th e 
case of the known Hot Jupiters fe.g. lMavor fc Udrvll2000h . 
However, new developments based on the estimate of the 
exospheric temperatur e of HP 209458 b, evaluat ed to be 
much higher than T e g IjVidal-Madiar et al.ll2003|) . call for 
a complete reconsideration of the question. 

These approaches are very interesting to explain the 
period distribution of planets very close to the central 
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stars (a<0.1AU, with a pile up around 3 days), espe- 
cially considering that the peak at very short periods 
is almost entirely formed by the lightest mass planets 
(m 2 sin i < 0.75 M Jup ). However, they do not explain the 
extended range of periods (P < 100 d) in which no planets 
more massive than ~ 2 Mj up are found. This large interval 
supports the non-migration of massive companions rather 
than the disappearance of close-in migrating massive plan- 
ets. 

4. Summary and concluding remarks 

To summarise the main observational features pointed out 
in this paper and their main implications, we have that: 

1) No massive planets {mi sini > ~ 2 Mj up ) are found 
on short-period orbits (P<~100d) around single stars. 
This is not an observational bias as these candidates are 
the easiest ones to detect. 

2) The maximum mass of detected planets per period 
interval increases with distance to the central star (Fig.EJ). 
This is also a solid result as massive planets found at a 
given distance are easier to detect closer in. 

3) The above points 1) and 2) suggest that the migra- 
tion rate of planets decreases with increasing mass of the 
planetary companion. This result agrees wit h recent simu- 
lations of m i grating planets in v iscous disks l|Trilling et alJ 
Il998l l2002t iNelson et alll200C|) . We thus expect a large 
number of massive planets to be on long-period orbits and 
so to be still undetected because of the limited duration of 
the present surveys. A large number a lower-mass planets 
probably also exist on long-period orbits; they are how- 
ever more difficult to detect. These represent primary tar- 
gets for future h igher-precision surveys (like e.g. HARPS, 
Pepe et al . 2002) whereas the youngest among the formers 
are interesting targets for direct imaging of planetary-type 
objects. This result is also very important for on-going 
radial-velocity planet searches from which we now expect 
an increasing number of planet candidates as the survey 
durations will grow. The "older" the survey, the higher the 
expected detection rate. In this context, the first epoch 
measurement of the targets is an important parameter for 
planet detection. 

4) We observe a shortage of planets with periods in 
the 10-100 d range. This valley in the distribution is lo- 
cated just in-between the peak of light planets that have 
migrated inwards and were stopped close to the central 
star and the rise of the distribution due to the increasing 
number of detected planets with longer periods. 

5) Up to now, no planet candidates with very low 
masses (ms sini < 0.75 Mj up ) have been detected on orbits 
with periods longer than ~ 100 days. This trend seems sig- 
nificant and not due to small-number statistics or to the 
particular shape of the period or separation distributions. 
The limit is sharp in mass, indicating a strong transition 
in the migration process for different mass regimes. These 
features are in agreement with predictions of runaway 
migration simulations of planets lighter than ~ 1 Mj up 



in m assive protoplanetary disks IjMasset &: Papaloizoul 
l2003h . 

Note however that this lack of light planets on inter- 
mediate periods may not be confirmed for much lighter 
planets or for planets in the considered mass range but 
formed far out from the central star (e.g. Neptune- and 
Uranus- type planets), that are still out of reach of the 
present observational facilities and detection techniques. 

In conclusion, we have pointed out a clear dependence 
of migration on mass and distance from the central star, 
providing important constraints for models of planetary 
formation (migration). More quantitative constraints for 
the different mass regimes will require a better detection 
capability of the surveys, to diminish the effect of the ob- 
servational bias in the distributions and thus increase the 
confidence we can put in the derived trends, especially for 
the lower-mass planets. 
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